Generalized Fluid System Simulation Program (GFSSP) 
Version 6 - General Purpose Thermo-Fluid Network 

Analysis Software 



Alok Majumdar, Andre Leclair, Ric Moore 
NASA/Marshall Space Flight Center 

& 

Paul Schallhorn 
NASA/Kennedy Space Center 

Thermal Fluids Analysis Workshop (TFAWS) 
August 15-19, 2011, Newport News, VA 


Content 



• Introduction 

• Additional Capabilities of Version 6 

— Fluid Mixture Option with Phase Change 

— Pressure Regulator Model with Forward Looking 
Algorithm 

— Prescribed Flow Option 
— Two-dimensional Navier-Stokes Solver 
- SI Option 

• Concluding Remarks 


2 


Introduction 
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• GFSSP stands for Generalized Fluid System Simulation 
Program 

• It is a general-purpose computer program to compute 
pressure, temperature and flow distribution in a flow 
network 

• It was primarily developed to analyze 

- Internal Flow Analysis of a Turbopump 

- Transient Flow Analysis of a Propulsion System 

• GFSSP development started in 1994 with an objective to 
provide a generalized and easy to use flow analysis tool 
for thermo-fluid systems 
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Development History 
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• Version 1 .4 (Steady State) was released in 1996 

• Version 2.01 (Thermodynamic Transient) was released 
in 1998 

• Version 3.0 (User Subroutine) was released in 1999 

• Graphical User Interface, VTASC was developed in 2000 

• Selected for NASA Software of the Year Award in 2001 

• Version 4.0 (Fluid Transient and post-processing 
capability) is released in 2003 

• Version 5 (Conjugate Heat Transfer) is released in 2007 
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Network Definition 
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Program Structure 
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Graphical User 


Solver & Property 
Module 


User Subroutines 
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Mathematical Closure 



Unknown Variables 

1. Pressure 

2. Flowrate 

3. Fluid Temperature 

4. Solid Temperature 

5. Specie Concentrations 

6. Mass 


Available Equations to Solve 

1. Mass Conservation Equation 

2. Momentum Conservation Equation 

3. Energy Conservation Equation of Fluid 

4. Energy Conservation Equation of Solid 

5. Conservation Equations for Mass Fraction of Species 

6. Thermodynamic Equation of State 
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Graphical User Interface 
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Capabilities 



• Steady or unsteady flow 

• Compressible or incompressible flow 

• Single fluid or mixture 

• 25 flow resistance and 33 fluid options 

• Options for new components and physics through User 
Subroutine 

• Options for new fluid through table look-up 

• Conjugate Heat Transfer 

• Interface with Thermal Analysis Code, SINDA-G/PATRAN 

• Translator of SINDA/Fluint Model 
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Additional Capabilities of Version 6 



Fluid Mixture Option with Phase Change 

Pressure Regulator Model with Forward Looking 
Algorithm 

Prescribed Flow Option 
Two-dimensional Navier-Stokes Solver 
SI Option 



Fluid Mixture Option with Phase Chang 



• The mixture capability in earlier versions of GFSSP does not 
allow phase change in any constituent of the mixture 

• In liquid propulsion applications, there are situations where 
one of the constituents is saturated, i.e. mixture of liquid and 
vapor in equilibrium 

— For example during purging of liquid oxygen by ambient helium, a 
mixture of helium, L02 and G02 exist 

• Why is there such a limitation? 

— Because the energy conservation equation of the mixture is solved in 
terms of temperature 

— For calculating phase change, energy equation for each species must 
be solved in terms of enthalpy or entropy 
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Mathematical Formulation 



• Mass Conservation 

- Mixture Mass 

- Concentration of Species 

• Momentum Conservation 

- Mixture Momentum 

• Energy Conservation 

- Temperature option 

• Energy Conservation is formulated in terms of temperature 

• Applicable for gas mixture 

- Enthalpy option - 1 

• Temperature is calculated by an iterative Newton-Raphson method 

- Enthalpy option - 2 

• Separate Energy Equations are solved for Individual Species 

• Applicable for liquid-gas mixture with phase change 
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Mixture Enthalpy Equation 
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Temperature Equation 
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Temperature equation is solved iteratively adjusting T, 
until right hand side of Temperature equation becomes 
zero 
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Thermodynamic Properties 



• Temperature and other properties of individual species are calculated from 


node pressure and enthalpy of the species: 


• The nodal properties are calculated by 
averaging the properties of species: 


n r _ 

Pi = Pik 

k = 1 


n f _ 

Pi = Yj Cik Pik 

k = 1 


T ik =f{Pi’ h n<) 
=f{PiAk) 
Pik =f{PiAk ) 
K ik =f(Pi’ h ik) 
C P ,„ =f(Pi> l >ik) 


• Temperature is currently 
calculated by averaging based on 
molar concentration of species 

• Alternate method of 
temperature calculation based 
on Vapor Liquid Equilibrium for 
multi-component, multi-phase 
mixture is in progress 
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Pogo Accumulator with Charge Line and Dump line 
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GFSSP Model of Pogo Accumulator & Drain Line 
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Charging of helium & draining of He-L02 mixture 





Displacement of Oxygen from 
the Accumulator during 
charging 





Pressure Regulator Model with Forward 

Looking Algorithm 



In Marching Algorithm, area is 
guessed and adjusted only once 
in each time step 

Adjustment of area is calculated 
based on difference between 
calculated and desired pressure 

Area adjustment can be done by 
backward differencing algorithm 
(Schallhorn-Majumdar) or 
forward looking algorithm 
(Schallhorn-Hass) 

Schallhorn-Hass Algorithm has 
been implemented in GFSSP 
Version 602 
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Backward & Forward Differencing Algorithm 
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Application of Forward Looking Algorithm 
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Reference: Forward Looking Pressure Regulator Algorithm for Improved 
Modeling Performance with the Generalized Fluid System Simulation 
Program by Paul Schallhorn & Neal Hass, AIAA Paper No. 2004-3667 
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Fixed Flow Option 
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• A new branch option has been introduced to fix flowrate in a 
given branch 

• The fixed flow branch can only be located adjacent to a 
Boundary Node 

• For unsteady option, a history file will be needed to specify 
flowrate and area at all timesteps 

• With this new option a user can prescribe either pressure or 
flowrate as boundary condition 

• Flow Regulator option is also available in unsteady mode to fix 
flowrate in an internal branch 
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Algorithm for Fixed Flow Option(Schallhorn) 
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New Resistance Option - Fixed Flow 
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Properties of Fixed Flow Option 
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Results of Fixed Flow Option 

Assigned Flowrates 
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Two-dimensional Navier-Stokes Solver 


Higher fidelity solutions are 
often needed that are not 
within the capacity of system 
level codes. 

GFSSP's momentum equation 
has been extended to perform 
multi-dimensional calculation 


Top Wall Top Wall Top Wall Top Wall Top Wall Top Wall 


► ► ► ► ► 



7 

CO 

14 

co 

21 

CO 

28 

CO 

35 

CO 

42 

CO 

49 

Side Wall 

c4> 




s 


c4> 




c4> 


be 

(4849) 


6 

co 

13 

co 

20 

CO 

27 

CO 

34 

CO 

41 

CO 

48 

Side Wall 



I 


(j92() 


bb 

(2627) 




be 

(4041) 


be 

(4748) 


5 

co 

12 

co 

19 

CO 

26 

CO 

33 

CO 

40 

CO 

47 

Side Wall 

c«> 






(2526) 


<s> 


be 

(3940) 


Cs> 


4 

CO 

11 

co 

18 

CO 

25 

CO 

32 

CO 

39 

CO 

46 

Side Wall 





(j71j) 


<^> 


r^i 


be 

(3839) 


be 

(4546) 


3 

-co 

10 

O' 

17 

CO 

24 

CO 

31 

CO 

38 

CO 

45 

Side Wall 





s 


<^D> 


c^> 


c§> 


be 

(4445) 


2 

-co 

9 

CO 

16 

CO 

23 

CO 

30 

CO 

37 

CO 

44 

Side Wall 





s 




(293t) 


<^5 


be 

(4344) 


1 

-CO 

8 

CO 

15 

CO 

22 

CO 

29 

CO 

36 

CO 

43 


Bottom Wall Bottom Wall Bottom Wall Bottom Wall Bottom Wall Bottom Wall 


Side Wall 


Side Wall 


Side Wall 


Side Wall 


Side Wall 


Side Wall 



< 

12 inches 


Shear Driven Square Cavity Centerline Velocity Distribution 



Velocity Field and Pressure 
Contours 



Dimensionless Velocity 
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NodeProperties -- SI UNITS 
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• SI Option is for 
input/output 

• GFSSP solver 
works in Engineering 
Unit 

• User Subroutine 
must be in 
Engineering Unit 
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Concluding Remarks 
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• GFSSP Version 6 will have additional capabilities to model: 

- Fluid Mixture Option with Phase Change 

- Pressure Regulator Model with Forward Looking Algorithm 

- Prescribed Flow Option 

- Two-dimensional Navier-Stokes Solver 

- SI Option 

• GFSSP is available (with no cost) to all Federal Government 
Organizations and their Contractors 

• Concepts/NREC has the license for commercial distribution to 
domestic and international Companies or Universities 

• A process is in work to make an educational version available to all 
Accredited US Universities for teaching and research 
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